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Abstract 

In this paper we study the cumulative conditional expectation function (CCEF) in the copula con¬ 
text. It is shown how to compute CCEF in terms of the cumulative copula function, this natural 
representation allows to deduce some useful properties, for instance with applications to convex com¬ 
bination of copulas. We introduce approximations of CCEF based on Bernstein polynomial copulas. 
We introduce estimators for CCEF, which were constructed through Bernstein polynomial estimators 
for copulas. The estimators are asymptotically normal and biased for CCEF. 

1 Introduction 

In this paper we explore E[U|t/ < m] as a function of it G (0,1), we denote this quantity as cumulative 
conditional expectation function. One of the motivations of working with this quantity is its use for 
making decisions in real problems. Our target is to give the foundation for the construction of approx¬ 
imations and estimators for the cumulative conditional expected function. Also, we present in the first 
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case the rate of approximation and in the second, the asymptotic distribution of the empirical process 
related to the estimator. 


2 Cumulative conditional expectation 


For some applied problems, it is helpful knowing the mean performance of one variable conditioned 
to another variable restricted to a region and not just to a single observation. In order to reach that 
knowledge we study the following measure. 

Definition 1 Let {U, V) he a random vector with associated 2-copula C. Then the cumulative condi¬ 
tional expectation function ofV on U is given by Rc{u) = E[y|{7 < u], Vm S (0,1). 


Remark 1 The cases u = 1 and u = 0 are excluded because they restrict themselves to the copula 
regression function. Indeed E[V|17 < 1] = E[V] = 0.5 is trivial and E[V|[/ < 0] = E[V|[/ = 0] has been 
studied in m- 

We note the previous concept was already explored in [3, but that work focused on some bounds of 
the measure Re- The following theorem gives a representation for Rc{u) that let us to simplify its 
computation. 

Theorem 1 Let {U, V) he a random vector with associated 2-copula C. Then Rc{u) = 1—^ fg C(u, v)dv, 
VmG (0,1). 


In order to simplify the notation, hereafter we denote the partial derivatives as done in [^, this is 
C(2)(n,n) = and It can be 


seen that E[V|[/ < u] is the average of the copula regression function E[I/|{7 = u]. 


Corollary 1 Under the hypotheses of theorem{^ Rc{u) = ^ /q“E[V|[/ = w\dw, Vu G (0,1). 
Proof. By Theorem[T]we have Rc{u) = ^ ^ Jq C{u,v)dv, then 


Rciu) 




w, v)dwdv 


— J (l — J C^^\w,v)dv'jdw. 
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The last term is equal to ^ 'E\V\U = w\dw according to [13]. ■ 

By means of Theorem[T]we can compute the cumulative conditional expectation function for several 
cases, as it was explored in |5|. The following result shows that the cumulative conditional expectation 
of a convex mixture of copulas is the convex mixture of the cumulative conditional expectations. 

Corollary 2 Under the hypotheses of theorem]^ let C be the convex eombination copula, C(u,v) = 
where is a copula function, 0 < p.y < 1 for 1 < 7 < T and = 1- Then, 

Rc{u) = E!y=iP 7 i?c.,(u). 

In the next result we study a specific class of copulas, the polynomial ones. Recalling that a copula 
has polynomial cross sections in u if it can be written as C{u, v) = oii{v)u'‘ for each u S [0,1], to 
suitable functions ai,i = 1,..., fc. More details about these copulas can be found in [10] §3.2. As an 
immediate consequence of Theorem [I] we have the following corollary. 

Corollary 3 Under the hypotheses of theorem\^ if C{u,v) = Oii{v)u’'. Then, Rc{u) is a poly¬ 

nomial of degree k — 1 given by Rc{u) = 1 — Yli=o fo 

Example 1 Quadratic and cubic cross sections copulas. 

1. The Farlie-Gumbel-Morgenstern (FGM) family contains all the copulas with quadratic cross sec¬ 
tions in both variables. C(u,v) = uv -\- 6uv(\ — u)(l — v) with 9 € [—1,1], then Rc(u) = 
(3 + 6l(u- l))/ 6 . 

2. The Lin’s iterated FGM family is defined by C{u, v) = uv-\- 6uv{l — u)(l — v) (l + ip{l — it)(l — u)) 

with 9 andifi real numbers satisfying 9 G [—1,1] and—1—9 < 9{l-\-Lp) < (3—6*+(9— 6 d—3d^)^/^)/2. 
C has cubic cross sections in both variables. For u € [0,1], C{u,v) = ctiiv)u'‘ where 

ai{v) = {1-\-9-\-(p9)v-\-{—9 — 2ip9)v'^-\-g)9v'^ 
a2(v) = (—9 — 2<p9)v-\-(9-\-4(p9)v'^ — 2ip9v^ 
o^siv) = ‘p9v — 2ip9v^ + q:}9v^. 

Then, Rc{u) = (i - | 
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3 Bernstein Copulas 

Definition 2 Given a 2-copula C and a grid of points (uk, vi) € [0,1]^ with k,l = the 

Bernstein copula approximation of C of order m is given by 

mm j j 

BmC{u,v) ='^'^c(—,—)Pk,m{u)PL,n{v), V(lt, t;) G [ 0 , 1 ]^, 

\m mJ 


where Pj^m{x) = (™)a;-^(l — a;)™ ^, x G [0,1]. 

The function B^C approximates the copula C and is a copula itself. The usefulness of the next 
theorem is to numerically simplify the computation of Rc, through solving polynomial integrals, when 
C is a parametric copula described by a complex analytical expression. 

Definition 3 Consider a 2-copula C, and a grid of points (uk, vi) G [0,1]^ with k,l = Set 

the Bernstein approximation of Rc of order m as 


RBrr,c{u) = 1 — 


1 


(m + l)u 


mm J J 

\m mJ 




where Pk,m{x) = {X)x^i^ - x)”^ x G [0,1]. 

Theorem 2 Under the hypotheses of theorem]^ Rb^c{u) converges to Rc{u) pointwise, when m —>■ 
oo. 


For appropriate copulas it is possible to give a rate of the approximation for Rc, and also we observe 
that this rate is improved when u increases. 

Theorem 3 Let C be a continuous copula with first order partial derivatives being Lipschitz and Uq G 
(0,1). For u G [mq, 1) it holds \Rc{u) — Rb^c{u)\ < for a constant M. Then, Rb^C converges 

uniformly to Rc in [uq, 1), when m —>■ oo. 
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4 Estimation 


Consider a bivariate random sample of the vector {U,V) with associated unknown 2- 

copula C. Let be Ru^ = rank of C/j in {t/i,..., and = rank of 1^- in {Vi,..., 14}. Denote the 
empirical copula as 

Cn{u,v) = - ^ < u)l(—^ <v), Vm,W e [0, 1]. 

i=i 


By using the Bernstein estimator of the copula function C proposed in m we obtain the estimator of 
the cumulative conditional expectation function. 

Definition 4 Given a 2-copula C, the estimator of Rc{u) is 


<?„(«) = 1 — 


Vug (0,1), 

\m m/ 


[m - 


l)u 


k=0 1=0 


where C„ is the empirical copula. 

In order to explore the asymptotic behavior of the estimator RB^Cn (^) make use of some known 
results about the Bernstein estimator. In [3] (Theorem 2 and Remark 3) it is proved that for a copula 
C with bounded third order partial derivatives on (0,1)^, if y/nim —)> d, 0 < d < oo, when n —>■ oo, 
then the process •/ri^i3mC'„(u, u) — (7(u, ^ Qc{u,v) in the space Z°°((0,1)^) of bounded functions, 

where ^ denotes weak convergence. In addition, the limiting process Qc{u,v) is a tight Gaussian 
process with mean function dh(u, v) where 


b{u,v) = i ^u(l — u)C^^'^\u,v) + u(l — u)C^^’^^(u, 


( 1 ) 


and covariance function given by E\h{u,v)h{u',v')] for 0 < u, u', v, v' < 1 with 


h{u, v) = I{U <u,V <v)- C(u, v) - C^^'> {u, v) {l{U <u)-u)- (u, u) {l{V <v)-v). (2) 
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This time we consider the process i fg C/c(u, v)dv, with u G (0,1), which is the process that defines 
the asymptotic behavior of the proposed estimator of Rc- 

Theorem 4 Let C be a 2-copula with hounded third order partial derivatives on (0,1)^. If n and m 
are positive integers such that y/nlm —>■ d with 0 < d < co when n,m ^ oo, then for u G (0,1), 


RsmcAu) - Rc{u )]— / Qc{u,v)dv. 

> u Jo 

~u fo Gciu,v)dv is a Gaussian process with mean function d{^ — Rc{u)) + fg C^^^^\u,v)dv 

and variance function 


Var Qc{u,v)dv = -4 :Rq{u)[4 --)Rc{u) - 4(2 --)Rc{u)Ri{u,u) 

u Jo J u u 

+2(3-m) — 4(l- )'HI{u,u) — 2-\ -1 —- f {C^‘^\u,v))^vdv. 

XI XL XI Xi^ Jq 

where 2'Hi(u,m) = fg C^^^(u,v)dv. 


Remark 2 The covariance function of the process fg Qc{u,v)dv is 


Cov 


Qc(u,v)dv, - - / Qc{u',v')dv 


= —- f [ C(u/\u',v/\v')dv'dv-\-(l — Rc{uf){^Rc(u') — l')-\-TLo{u,u')+TLo{u',u) 

UU Jo Jo 

+ ( —7 ~ l) f C^^\u',v)dv [ C^'"\u,v)dv — R2 {u,u') — R2{u', u) 

u' 'Jo Jo 

H- [ [ C^'^\u' ,v')C^‘^\u,v)(v /\v')dvdv' — Rc{u)Rc{u') -\-Ri{u,u') -\-Rx{u' ,u). 

uu' Jo Jo 
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with 


'Hi{wi,W2) 

'H2iwi,W2) 

'Hz{wi,W2) 


1 


W1W2 Jo 


C'^^\wi,v)dv / C(wi,v)C‘'‘^\w2,v)dv 


—— [ [ C^^\w2,v)C{wi,v Av')dvdv' — (1 — Rc{wi))Rci W2) 
W1W2 Jo Jo 

1 

( --- (Rc(wi A W2) - 1) + 1 - 2 Rc(wi)) / C^^^^(w2,v)dv. 

\W1VW2' j Jo 


for wiV W 2 = max{'u;i, ^ 2 } and wi Aw 2 = min{wi, ■u; 2 }- 
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